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Abstract: We propose a QCD motivated theoretical approach to high energy soft interactions, which 
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that the survival probability for the diffractive Higgs production at the LHC energy is small (less than 
1 1%), and investigate the influence of suggested corrections e.g. threshhold effects and semi-enhanced 
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1. Introduction 

For the past three decades the physics of soft interactions has been considered an unsavory topic. Indeed, 
the relation between high energy phenomenology and the microscopic theory: QCD, has been neglected, 
as there has been no progress in solving the problem of the confinement of quarks and gluons. The goal of 
this letter is to present an attempt to construct a self consistent theoretical approach to strong interactions 
at high energy, which provides a natural bridge to the perturbative QCD description at short distances. 

Following on the results of our paper [1], we question the widely held prejudice that high energy 
scattering stems from long distances. We believe, that the only dimensional scale that governs the high 
energy interaction, is the saturation momentum (Q s ), which appears both in pQCD at high parton density 
[2,3,5,6], and in the non-perturbative approach to QCD based on the Colour Glass Condensate [4,7]). 
The saturation momentum increases at high energies and, therefore, we believe that the higher the energy, 
the shorter are the distances involved in the interaction. The key argument supporting our hypothesis is 
the fact that in fits to data, the slope of the Pomeron trajectory {a'p) is very small [1,8,9]. At first sight 
a small a'jp contradicts one of the principle results of Regge theory, i.e. that all resonances lie on Regge 
trajectories. However, the first theoretical analysis of the Pomeron's structure in a theory in which we can 
treat the region of strong couplings (the strongly coupled N=4 super Yang-Mills theory with AdS/CFT 
correspondence), given in Ref. [10], shows that a'jp — > for t < 0. This result is compatible with our 
approach. In the same theory a'jp > for t > 0, which is in reasonable agreement with Regge theory and 
high energy phenomenology (see Fig. [l] which was taken from Ref. [10]). Therefore, the small value of a'jp 
in the scattering region is not only possible, but is a requisite, so as to be consistent with a strong coupling 
theory. 
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In perturbative QCD a'jp for a BFKL Pomeron 




[11], which determines the high energy scatter- 
ing amplitude, is equal to zero, since a'jp oc l/Q 2 ■ 
0. W = ^/s is the energy in the cm. frame. 
Therefore, the assumption that a'jp — > 0, provides 
the necessary condition so as to link the strong 
interactions with the hard interactions that are 
described by pQCD. 
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The only microscopic explanation of Regge 
theory is provided by the partonic approach [13, 
14], in which the slope of the Pomeron trajec- 
tory is related to the mean transverse momentum 
of the exchanged partons, a'jp oc 1/ < p t >. 
< pt > is the mean parton momentum. The corn- 



Figure 1: The behaviour of the Pomeron trajectory in N=4 monly held view from high energy phenomenol- 
SYM accordingly to Ref. [10]. The figure is taken from Ref. ogy, is that a'jp = 0.25 GeV~ 2 [15]. 



are small, and the shrinkage of the forward elastic peak is induced by the slope of the Pomeron trajectory. 
The above supposition is not compatible with the high value of ajp(0) implied by the ISR-Tevatron cross 
section data. In Ref. [1] an alternative option was suggested in which eikonal shadowing corrections are sig- 
nificant yielding a diminishing a'jp = 0.01 GeV~ 2 value. This result is not affected by a more complicated 
form of the Pomeron interactions. A similar result was obtained, also, in Ref. [8] in which the experi- 
mental elastic data was reproduced with a'jp = *. We conclude that a small a'jp is a reasonable, data 
compatible, option which induces a large typical parton momentum. The resulting running QCD coupling 
a s = 7r/Mn ^< pi > /-^Qcj^j <C 1. We shall consider it as the small parameter in the pQCD estimates of 
the Pomeron-Pomeron interaction vertices. We wish to stress that, regardless of its intuitive appeal, the 
parton model does not provide a reliable approach, since it presumes a short range interaction in 
rapidity, while the exchange of gluons in QCD generate long range rapidity interactions. Nev- 
ertheless, the model convincingly illustrates, that the smallness of a'jp, implies that the typical interaction 
distances are short. As we have discussed, we estimate that they are of the order of 1/Q S . 

The second assumption underlying our approach is that only the triple Pomeron interaction is essential 
at high energy. This does not precludes multi Pomeron vertices which are estimated by fan diagrams. As 
such these vertices are not local in rapidity. This has been proven in pQCD (see Refs. [17,18]). This selection 
rule (considering only the triple Pomeron vertex) postulates that in the leading log x approximation of 
pQCD, with a large number of colours (N c ^> 1), all other possible interactions are suppressed. It is well 
known that at high energy, and for N c ^> 1, the correct degrees of freedom are the colourless dipoles [19]. In 
terms of these dipoles, the selection rule reads, that only the decay of one dipole to two dipoles contributes 
to the high energy amplitude. This implies that a colourless dipole does not interact with the target via 
*In their last five papers (see Refs. [9]) the Durham group has considered a diminishing, non zero, a'p. 
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The phenomenological value of a'jp stems from 
the supposition, that the shadowing corrections 
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Figure 2: The different contributions to Pomeron Green's function ( enhanced diagrams of Fig. ||-a) and the semi- 
enhanced diagrams that has not been taken into account in the Pomeron Green's function. 

multi Glauber-type rescatterings, but only through the splitting of one dipole into two dipoles and the 
merging of two dipoles to one dipole. This postulate provides a natural bridge between our theoretical 
approach to high energy scattering, and the corresponding pQCD high energy scattering amplitude. 

In this paper we present a theoretical approach for soft interactions at high energies, incorporating the 
following assumptions: 

1. a'jp = 0. 

2. The Pomeron interactions reduce via fan diagrams to triple Pomeron vertices. 

3. a'jft ^ and it's slope, as well as other parameters of the Reggeons, should be determined from 
experiment. 

4. A Pomeron - hadron interaction is treated phenomenologically using the eikonal formulae, in the 
framework of a two channel formalism, based on the Good- Walker mechanism [20]. 

The first two assumptions allow us to build a theory of interacting Pomerons, while the last two points 
are phenomenological. Details of the two channel (Good- Walker mechanism, sec Rcf. [20]) that we use for 
describing the elastic and diffractive scattering eigenstates are given in Ref. [1]. 

2. The Pomeron interaction (Pomeron Green's function) 

Our key idea is to use the general results of pQCD relating to high energy scattering [12]. Consequently: 
1) In the leading order approximation of pQCD, only one Pomeron splitting into two Pomerons, and two 
Pomerons merging into one Pomeron are taken into account [17, 18]. All other Pomeron vertices do not 
explicitly appear in the leading log x approximation of pQCD. We, therefore, restrict ourselves to sum only 
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Pomeron diagrams containing triple Pomeron vertices. 

2) We consider a theory with a'jp = 0. 

3) Using our approach, we can predict the anticipated values of all Pomeron interactions: the intercept of 
the Pomeron above unity Ajp is oc as, and the triple Pomeron vertex coupling g^jp oc a|. 

The theory which includes all the above ingredients can be formulated in terms of a generating function 
[19,26] 

Z(y, u) = V P n (y) u n , (2.1) 



n 



where, P n (y) is the probability to find n-Pomerons (dipoles) at rapidity y. At rapidity y = Y = ln(s/so) 
we can impose an arbitrary initial condition. For example, demanding that there is only one fastest parton 
(dipole), which is P\{y = Y) = 1, while P n >i(y = Y) = 0. In this case we have the following initial 
condition for the generating function 

Z(y = Y)=u. (2.2) 

At u = 1 

Z(y, u = 1) = 1, (2.3) 

which follows from the physical meaning of P n as a probability. The solution, with these two conditions, 
will give us the sum of enhanced diagrams (see Fig. §-a). 

For the function Z (u) the following simple equation can be written (see Ref. [1] and references therein) 

dZ(y, u) dZjy, u) d 2 Z(y, u) 
d^T~ = -r(l^2)u(l - u) Qu + r(2->l)u(l - u) Q2u , (2.4) 

where, T(l —* 2) describes the decay of one Pomeron (dipole) into two Pomerons (dipoles), while T(2 — > 1) 
relates to the merging of two Pomerons (dipoles) into one Pomeron (dipole). 

The first term on the r.h.s. describes two processes: the term proportional to u can be viewed as a 
probability of a dipole annihilation in the rapidity range (y to y — dy) (death term). The term with factor 
u 2 is a probability to create one extra dipole (birth term). The second term on the r.h.s. describes the 
same type of processes for Pomeron (dipole) merging. 

The description of the parton system given by Eq. (|2.4| ), is equivalent to the path integral formulation 
of the Pomeron interaction (see Ref. [1]). The path integral formulation of the Pomeron interaction leads 
to the automatic inclusion of i-channel unitarity constraints, while in the generating function formulation, 
s-channel unitarity has been taken into account. 

Using the functional Z, we find the scattering amplitude [6,26], using the following formula: 

N(Y) lraA d {Y) = f^t^l d l^l u) \ u=lln (Y = Y ,b), (2.5) 

n=l 

where, "y n (Y = Yoify 1S the scattering amplitude of n-partons (dipoles) at low energy. These amplitudes 
depend on the impact parameters which are the same for all n partons since a'jp = 0. Consequently, 
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we neglect the diffusion of partons in impact parameter space. Eq. (|2.5| ) corresponds to the partonic 
approach [13,14], in which a high energy scattering process can be viewed as a two stage sequence. The first 
stage is the development of the partonic wave function, which we consider by introducing the generating 
function Z . The second stage is given by Eq. ( |2.5| ) corresponds to the interaction of the lowest energy 
partons ('wee' partons) with the target, which is described by the amplitudes j n (Y = Yo,b). Assuming 
that there are no correlations between the interacting partons (dipoles) at low energy, we take ~f n (Y = 
Y ,b) = 7 ?(y = Y ,b) [6]. 

The generating function approach given by Eq. fl2.1|), Eq. fl2.4| ) and Eq. (|2.5|), has the advantage that 
it can be solved analytically (see Ref. [29]). This solution leads to a constant cross section at high energy, 
while the interaction without the four Pomeron term, decreases at high energy [27] (see also Refs. [25,28]). 
The problem is that the inclusion of the four Pomeron term contradicts our QCD motivated approach. 
This highlights the importance of s-channel unitarity studies of the asymptotic high energy behaviour of 
the scattering amplitude. The above problem is not relevant to our present investigation as we restrict 
the energy range where we can trust our approach, to energies at which the four Pomeron term is not yet 
significant. As discussed in our paper [1] this range is W < 10 5 GeV. 

We wish to emphasis that the form of the equation, as well as the functional formulation of the Pomeron 
interaction, is the same as that for interacting dipoles in QCD. The only difference is that in QCD the size 
of the dipole can be changed due to an interaction, while in our case the size of the dipole is fixed. The idea 
that the transition from short distances to long distances goes through the stage of freezing the typical size 
of the interacting dipoles, has been around for a long time, and has both theoretical and phenomenological 
justification (see, for example, Refs. [21-24]). This gives us further confidence that our approach provides 
a natural link with the perturbative QCD approach. 

Unfortunately, the exact solution to Eq. (2.4) with arbitrary initial and boundary conditions, is not 
very transparent. In a number of papers [1,30,32-34,37], it has been shown that, using Mueller-Patel- 
Salam-Iancu approximation, one can simplify the solution, and reduce it to a simple analytical form. This 
approximation yields a solution with an accuracy of the order of 7 exp (— Ajp Y) oc a| exp (—asY) 1. 

In ref. [1], and in this letter, we only calculate the Green's function of the Pomeron, neglecting for 
example, diagrams of the Fig. g-b type. The reason for this procedure is that the interaction of the Pomeron 
described by the exact Green's function, will be simpler than the interaction of the bare Pomerons (see 
Ref. [35] in which this approach was suggested). 

The exact expression for the Pomeron Green's function has the form 

G F (Y) = l-exp^^^O,^), (2.6) 

where T (0,x) is the incomplete gamma function (see 8.350 - 8.359 in Ref. [38]) and 

T{Y) = 7 e A ^ y . (2.7) 

where, 7 has a simple meaning of being the amplitude of of the two dipoles interaction at low energy. The 
derivation of Eq. (2J3) as well as derivation of the formulae for the cross sections of the single and double 
diffraction production, are given in Ref. [1]. 
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3. The results of the fit 



The pertinent details of our fit to the experimental data, and our determination of the relevant parameters 
of the model, needed to describe the soft interactions, are contained in [1]. In this section we will only 
mention the salient features, and results of the fit. 



o tot (s)(mb) 



K -I 




W(GeV) 

Figure 3: The total cross section [a tot = l/2(&tot(pp) + 
(Ttot(pP))) at low energies. The curve illustrates our 
parametrization. 



Our fit is based on 55 experimental data points, 
which includes the p-p and p-p total cross sections, 
integrated elastic cross sections, integrated single 
and double diffraction cross sections, and the for- 
ward slope of the elastic cross section in the ISR- 
Tevatron energy range. The model gives a good 
reproduction of the data, with a x 2 /d.o.f. ~ 1.25. 
A significant contribution to x 2 /d.o.f. stems from 
the uncertainty for the value of two single diffrac- 
tion cross sections, and of the total cross section 
at the Tevatron, without these three points we 
achieve a x 2 /d.o.f. = 1.0. In addition to the 
quantities contained in the data base, we obtain a 
good description of the CDF [39] differential elas- 
tic cross sections and the single diffractive mass 
distribution at t = 0.05 GeV 2 . An impor- 
tant advantage of our approach, is that the model 
provides a very good reproduction of the double 
diffractive (DD) data points. Other attempts to 
describe the DD data e.g. (see Refs. [9]), were not 
successful in reproducing the DD experimental re- 



sults over the whole energy range. 

The predicted energy behaviour of the above mentioned cross sections, as well as the forward slope of 
the elastic cross section, are given in detail in Ref . [1] . We mention here only our results for the LHC energy, 
which are of particular interest i.e. at W=14 TeV we have a tot = 92.1 mb, a e \ = 20.9m6, a so \ = 11.8m6, 
a dd = 6.08 mb and B el = 20.6GeV~ 2 . 

We note that the lowest energy data point in our fit [1] to the total cross section was at W = 19.2 
GeV. In Fig. 3, we show the results of our calculations in which we have extrapolated our model to lower 
energies. As can be seen, it provides a fair description of the data even at energies as low as W = 7 GeV, 
where we underestimate the experimental numbers just by about 10%. 

4. Survival probability for diffractive production of the Higgs boson. 



In our approach, the survival probability for the diffractive production of the Higgs boson, can be written 
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Tevatron 


LHC 


W = 


itfGeV 




GLMM 


KMR 


GLMM 


KMR 


GLMM 


KMR 


(Ttot( nib ) 


73.29 


74.0 


92.1 


88.0 


108.0 


98.0 


<7 e ;(mb) 


16.3 


16.3 


20.9 


20.1 


24. 


22.9 




9.76 


10.9 


11.8 


13.3 


14.4 


15.7 


Jow M 

a sd 


8.56 


4.4 


10.52 


5.1 


12.2 


5.7 


high M 


1.2 


6.5 


1.28 


8.2 


2.2 


10.0 


o dd (mh) 


5.36 


7.2 


6.08 


13.4 


6.29 


17.3 


(a e i + a s d + add) jo-bat 


0.428 


0.464 


0.421 


0.531 


0.412 


0.57 


S 2ch(%) 


3.2 (2.7) 


2.7 - 4.8 


2.35 (0.45) 


1.2-3.2 


2.0 


0.9 - 2.5 


S enh(%) 


28.5 (100) 


100 


6.3 (100) 


100 


3.3 


100 


S 2 {%) 


0.91 (2.7) 


2.7 - 4.8 


0.15 (0.45) 


1.2-3.2 


0.066 


0.9 - 2.5 



Table 1: Comparison of the GLMM (this paper) and KMR [8] models. The numbers in parenthesis are the estimates 
for the contribution of the semi-enhanced diagrams to the survival probability. 



as a product of two factors, 

(| S 2 |) = (| S 2 nh |) (enhanced diagrams) x (| |) (Good-Walker mechanism) . (4-1) 



The second factor has been determined in a number of papers, and all the groups obtain approximately 
the same values. The first factor is new and was first evaluated in Ref. [1], where the exact formula for 
(| S 2 nh |) was derived. It has the form 

(\S 2 enh (MPSI)\)(Y) = (4.2) 

= S(T(Y,Y')) ~ 2e- A ^ Y - SY ^ 2 S 1 (T(Y,Y')) + e^^ 6 ^ 2 S 2 (T(Y,Y')) ; 
S(T) = ^{-T + (l + T) e^L^O,! 

Si(T) = ± j-r(l + T) + (l + 2T)e£r(o,^) } : (i.:>,) 



1 



S 2 (T) = ^ |T [(T - l) 2 - 2] + (1 + 3T) L [0, ~j j , (4.4) 
where, 

T(Y,Y') = 7 fe A -( y - y ') - l) (e^ Y ' - l) . (4.5) 



i 



1 



This factor diminishes the value of the survival probability, and causes a decrease in its magnitude by a 
factor of 5-6, in going from Tevatron to LHC energies (see Table |l|). 

After hearing our report on the above results at the "Diffraction'08" conference, Alexei Kaidalov posed 
two questions, (also see his talk at International Symposium on Multiparticle Dynamics 2008). The first 
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concerned the so called, threshold effect, which suggests that one can only use the one Pomeron exchange 
approximation at sufficiently large rapidities (Y > 2). In our approach there is no need for an additional 
cutoff in rapidity, since the scattering amplitude is described by the sum of the contributions of both 
Pomerons and Reggeons. 

In Fig. | we show the scattering cross section at rather low energies. Our parametrization provides a 
good description of the experimental data also at low energies. Since in our parametrization we include 
the contribution of the secondary Reggeons, there is no need for an additional cutoff. The second question 
concerns the semi-enhanced diagrams of the Fig. ||-b type. Indeed, in our fit since g± and 52 turn out to 
be large, we could anticipate a large contribution. A self consistent approach which includes all Pomeron 
diagrams will be published in a separate paper, which is presently being prepared [36]. In the present 
paper we only calculate the exclusive contribution of the diagrams of type Fig. |2|-lo , given by the following 
formula (the proof of this formula will be given in Ref. [36]) 




(I I) = ^T^nX I d 2 b{ 

1 + g l T(b';m l )T{ 1 ,A,y) + g k T(\b - b'^m^T^, A,y) 

g l 9k T(V; m) TQb- b'\; rm) T( 7 , A, y) \\ 

1 + 5i r(6';m i )T( 7 ,A,y) + g k T(\b - V\; rm) T( 7) A, y) J J ' 

In Eq. T(j,A,y) is given by Eq. (|?7|), and 



2 



(4.6) 



777? h 

r(6,m i )= -^ifiKt) (4.7) 

where K\ is the modified Bessel function of the second kind. The vertices masses rrij as well as the 
values of the parameters 7 and A and the form of the hard amplitude Ah were taken from Ref. [1] . 

The numerical estimates, in which the semi-enhanced contribution (Eqn. (4.6)) was exclusively in- 
cluded, are given in parenthesis in Table |l[ The semi-enhanced diagrams were neglected in [1]. Our 
present calculation, with only the semi-enhanced diagrams, serves only as a reminder that our calcula- 
tion in Ref. [1] is incomplete. We expect, though, S 2 (semi — enhanced) to considerably reduce (| 5| ch | 
) (Good- Walker mechanism), especially at the LHC energy. Estimates which include both semi-enhanced 
and enhanced diagrams will be published soon [36]. 

In spite of the fact that we do not need to take the threshold effect into account, we estimate its 
influence on the value of the survival probability. We assume that we can use Pomeron exchange only 
for Y > yo, in this case the first enhanced diagram shown in Fig. |i| has the following contribution (see a 
detailed discussion in Refs. [33,34]) 



,4 (Fig. |) = ^ e -2A(y- w ) + ( 7 e- AY ). (4. 



Therefore, in order to introduce the threshold correction we need to change 7 — ► 7 exp (Ayo), and multiply 
the Pomeron contribution by the step function O (Y — yo). The first observation means that in order to 
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maintain the quality of the experimental data description obtained without the threshold effect, we need to 
increase the value of 7 in our parametrization by a factor exp (Ayo)- We, then, need to take into account 
Q(Y — yo). Doing so we obtain the following values for (| S 2 nh |) at yo = 0,1.5,2.3: for the Tevatron 
energy we obtain values of 0.285, 0.7, 0.99 respectively, and for the LHC energy the corresponding values 
are 0.06, 0.12, 0.19 respectively. Our conclusion is that the threshold effect is not needed in our approach. 
Regardless, its inclusion does not change the fact that the value of (| S 2 nh |) decreases by a factor of 5-6 in 
going from the Tevatron to the LHC energy. It should be stressed that since Y — y\~yi — « l/A w 3, 
the theoretical accuracy associated with this calculation is low. 





y% - « l/A 



Figure 4: The first enhanced diagram. 



5. Conclusions 

Y-yiRt l/A 

In this letter we present a consistent theoretical 
approach that (i) provides a natural bridge both 
with a N=4 SYM Pomeron with a'p = 0, and 
with the pQCD approach; (ii) gives an analytical 
expression for summing the Pomeron interactions; 
(hi) leads to a very good (x 2 /d.o.f = 1) descrip- 
tion of all available data on strong interactions at 
high energy and (iv) predicts survival probability to be less that 1% at the LHC. 

Within the framework of our approach, we still need to solve the problems of 1) accounting for the 
semi-enhanced diagrams summed together with the enhanced ones (see Ref. [36]); 2) to develop a theoretical 
approach for the Pomeron-hadron interactions; and 3) to formulate the Pomeron-Reggeon interaction on 
a more sound theoretical basis. 

To illustrate our achievements and problems, we compare our approach with the work of the Durham 
group (KMR), whose results were presented in six papers recently published [8,9]. Both approaches consider 
a'jp as being small. 

In both programs the Pomeron interaction was taken into account. The difference between the two 
approaches is that KMR made an ad hoc "reasonable" assumption, that the multi-Pomeron vertices have 
the following form, for the transition of n Pomerons to m Pomerons 



q" 

Urn 



nm\ n+m ~ 2 g N /2 = nm A n + m ~ 3 g 3]P /2. 



(5.1) 



No theoretical arguments or theoretical models were offered in support of this assumption, which certainly 
contradicts the pQCD approach [17,18]. We view this assumption as problematic, since it eliminates the 
hope of building an approach that matches with pQCD. In paper 4 of Ref. [9], the system of two equations 
(Eq.(26) and Eq.(27) of this paper) has been suggested for the opacity and it is stated that they follow 
from the parton model. Since there is no proof of this claim, and while we are inclined to believe the 
authors, that the equations could indeed be proved, we should still bear in mind, that the parton approach 
(see Refs. [13, 14] is based on the assumption that there is only a short range rapidity interaction between 
partons, while we, due to exchange of gluons in QCD, have a long range rapidity interaction. 
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KMR numerically solved Eq.(5.1) using an iteration method. The accuracy of the solution and its 
stability were not discussed. Since our approach is a particular case of the KMR approach with only 
the triple Pomeron interaction, we are surprised that the interaction method works, since our Pomeron 
Green's function G F oc ^2 n {— l) n n\ (7 exp(Ay)) 71 . Such a series cannot not be reproduced by an iterative 
procedure. Unfortunately, KMR do not provide any details of their numerical procedure. 

For single diffraction in the region of high mass, KMR encounter the same problem as we have in this 
paper, since they do not include the full set of semi-enhanced diagrams (an example of a neglected diagram 
is the second diagram of Fig. Q-b). 

Although we consider the KMR approach to be much less theoretically reliable than ours, the numerical 
results regarding the elastic and total cross sections are very close. We consider this fact as confirmation 
that most properties of the high energy interaction are related to the Pomeron interaction, but not to the 
details of the interaction. 

We disagree with KMR, regarding the predicted value of the survival probability for central exclusive 
diffractive Higgs production. In Ref. [8] the effect of enhanced and semi-enhanced diagrams have been 
neglected, with some supportive arguments of why it is reasonable approximation, however, in the last 
five papers [9] some of the diagrams have been included. The main difference is that KMR consider the 
threshold effect as being very important, and by choosing Y > 2.3, they obtain a suppression factor 0.83 
for the Tevatron energy, and 0.67 for the LHC. 

We have explained why in our approach, it is not necessary to include the threshold effect separately. 
We have shown that the threshold effect does not change the fact, that the value of the survival probability 
decreases by a larger factor in going from the Tevatron to the LHC energy. We are of the opinion that, KMR 
should, first, systematically introduce the Reggeon contributions, and only after that discuss the threshold 
effect, both in the calculation of the survival probability and in the scattering amplitude, especially for the 
case of single diffraction production. 
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